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The status of lattice calculations for heavy quark systems is reviewed, focussing 
on weak matrix elements for leptonic and semi-leptonic decays of heavy mesons. 
After an assessment of the main systematic errors, results for the decay constants 
fo and fs, the B parameter describing B — B mixing, the Isgur-Wise function 
and the spectroscopy of heavy baryons are discussed. 



y^-^ 1 Introduction 

o\ 

Q\ The physics of hadrons containing heavy quarks plays an important role in 

1—1 the study of some of the CKM matrix elements, and therefore serves to test 

O | the consistency of the Standard Model. Lattice simulations of QCD provide 

D a non-perturbative treatment of hadronic processes, and are thus capable of 

dealing with large strong interaction effects in weak decay amplitudes. 

Lattice QCD replaces space-time by a four-dimensional hypercubic lattice 
of size L 3 ■ T. The sites are separated by the lattice spacing a, which acts as 
an UV cut-off. One problem encountered in current simulations is that typical 
values of a -1 lie in the range 2 — 3.5 GeV. Therefore one expects that dis- 
cretisation errors ( "lattice artefacts" ) will distort the results already for charm 
^ physics. Also, b quarks cannot be studied directly, since their mass is above 

the UV cut-off. One way of addressing this problem is to subtract the leading 
If} discretisation error by employing an 0(a)-improved lattice action 1 . Quantities 

0^ can then be computed safely around m c h ar m and extrapolated to the b quark 

mass. Alternatively, one can use the "static approximation" and perform the 
Q_i simulation at infinite heavy quark mass, using the leading term of an expan- 

sion of the heavy quark propagator in 1/ttiq. Furthermore, b quarks can be 
formulated using a non-relativistic expression for the QCD action. These three 
methods provide complementary information for b quark physics on the lattice. 

The other main systematic errors which affect our results include using the 
Quenched Approximation, in which the effects of quark loops are neglected. 
Furthermore, matrix elements computed on the lattice are related to their con- 
tinuum counterparts via finite renormalisation constants, due to explicit chiral 
symmetry breaking induced by the fermionic lattice action. The numerical 
values of these matching factors are subject to large uncertainties. Finally, 
lattice estimates of dimensionful quantities are affected by uncertainties in the 
lattice scale, which arise from the fact that different quantities used to estimate 
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Table 1: Lattice estimates for the pseudoscalar decay constants from different collaborations 



Collab. Jb [MeV] 



fp [MeV] 



fB./fE 



fojh 



MILC 3 
UKQCD 4 

PCW 5 
BLS 6 
ELC 7 



148(5)(14)(19) 

lou -6 -19 

180(50) 

187(10)(34)(15) 
205(40) 



180(4)(18)(16) 
185 + _t t 4 ? 
170(30) 
208(9)(35)(12) 
210(15) 



1.13(2)(9) 
1.22 t 4 3 
1.09(2)(5) 
1.11(6) 
1.08(6) 



1.09(1)(4) 
1.18 + _l 
1.09(2)(5) 
1.11(6) 
1.10(4) 



cT l [GeV] give different results. 

Most of the results discussed here come from a simulation by the UKQCD 
Collaboration, using a lattice of size 24 3 • 48 at j3 = 6/g 2 = 6.2 for which 
a -1 = 2.9(2) GeV. 2 For propagating quarks an 0(a)-improved fermion action 
was used. Results for b physics were obtained using propagating heavy quarks 
and also the static appoximation. 



2 Leptonic B decays and B° — B° mixing 

The decay constant of a heavy-light pseudoscalar meson, fp, can be extracted 
from the matrix elements of the lattice axial current via 

(0\A l : u (0)\P) ~M P fp/Z A (1) 

where Mp is the pseudoscalar mass, and Za = 0(1) is the matching factor 
relating the lattice matrix element to its continuum counterpart. In Fig. 1 wc 
show an example for the scaling behaviour of fpy/Mp with increasing heavy 
quark mass. According to HQET, this quantity behaves like a constant, but, as 
Fig. 1 shows, there are large deviations from this scaling law 4 . The estimates 
for f B , fo and the ratios f B ,/f B , /d s //d from various collaborations are 
shown in Tab. 1. A weighted average of the results in the table yields 

f B = 168 ± 30 MeV, f D = 196 ± 20 MeV (90% C.L.) (2) 

The phenomenologically interesting quantity for the study of CKM matrix 
elements is the combination f B \[B B , which requires knowledge of the B pa- 
rameter B B , defined via B B = a s (/x) _2 / /3 ° (B° | O l (ii) | £°)/f f B M|, where 
Ol is the AB = 2 four-fermion operator. In a recent study, UKQCD presented 
results for B B and f B using the static approximation 2 . The results for the B 
parameter are 

B Bd =0.99 ±* + _l B Bb =1.01 + _i 11 (3) 
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Figure 1: UKQCD's data for the quantity fpy'Mp q s (M p ) 2 / /3 o plotted versus 1/M 



P 



lattice units. The solid line denotes the extrapolation of the four points using propagating 
quarks. The point at infinite meson mass is obtained using the static approximation. 



D 



1 2 

1/M P 

where the first error is statistical and the second systematic. It should be 
stressed that there is a further 15-20% uncertainty coming from the matching 
procedure between the lattice and continuum matrix elements. It is reasonable, 
however, to expect that systematic effects cancel partially in the dimensionless 
ratios 2 

1.16 1^ (stat) 1 2 (syst) (4) 

1.37 1 g (stat) 1 4 (syst), (5) 

These results can be applied to the ratio of B° — B° mixing parameters x s /xd 
given by 

Xs_ = tb ± Vb ± _ M Bs J% B Bs \V ts \ 2 
x d r Bd r\ Bd M Bd P Bd B Bd \V t d\ 2 ' 

where r Bd a denote the Bd and B s meson lifetimes, and r) Bd a parametrise short- 
distance QCD corrections. Assuming r) Bd = i) Bs and using r Bd = 1.53(9) ps, 
t Bs = 1.54(14) ps, as well as our lattice estimate eq. (5), we find 

x s /x d = (1.38 ±0.17) \V ts \ 2 /\V td \ 2 . (7) 

In conjunction with the experimental value Xd = 0.76(6) 8 , the above ratios can 
also also be used to predict x s , provided | ^°| 2 is constrained using global fits 9 . 



ipVM; [a(M p )/«(M B )] 6/33 



0.10 



0.05 — , 



M B 
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For a fixed value of the B parameter Bk , x s is then obtained as a function of 
Jb^Mb- Choosing B K = 0.8, we find x s ~ 14 for f B VB~B ^ 170 Me V 2 . 

3 Semi-leptonic B decays and the Isgur-Wise function 

Processes like B — > Dlvi and B — > D*lvt can be parametrised in terms of six 
form factors, e.g. 



(D\V»\B) 
yj Mb Mo 



{vB+v D )[ 1 h + {w) + {v B -v D )[ 1 h-(Lo) (8) 



where u> = vb ■ «d denotes the product of 4- velocities if the B and D mesons. 
The Heavy Quark Symmetry (HQS) relates the six form factors to one universal 
form factor, called the Isgur-Wise function, which is normalised at zero 

recoil, £(u = 1) = 1. The form factor h + (to) is related to via 

/i+H = (i + /3 + M + 7+ H)eM (9) 

where (3+(co) parametrises radiative corrections between HQET and full QCD, 
and 7+(cj) denotes (unknown) corrections in X/thq. 

In our lattice study we concentrate on pseudoscalar (P) to pseudoscalar (f) 
transitions and extract the lattice form factor h + (ui) utt from the matrix ele- 
ment of the temporal component of the vector current. We determine the 
Isgur-Wise function through the ratio 

/» + M latt (i + /Mi)) „w v , im 

which, as is argued in 10 < 11 J is a good approximation, provided that the in- 
dependence of lattice artefacts and in 7+(w) is weak. 

In UKQCD's lattice simulation using propagating heavy quarks, one ob- 
serves a weak dependence of R+(u>) on the heavy quark mass. Furthermore, 
repeating the analysis for the axial form factor 1ia 1 (u>) relevant for B — > D*lvi 
decays, gives R + {lo) / Ra 1 (u) ~ l 11 . These observations support the universal- 
ity of £(w): one obtains the same estimate for the Isgur-Wise function over a 
range of different heavy quark masses (i.e. for different flavours), and also if 
the spin structure of the outgoing meson is changed. 

The lattice data for £(u>) are obtained for light quark masses around the 
strange quark mass and can be extrapolated to the chiral limit. The extrapo- 
lated Isgur-Wise function can then be parametrised in the form 

&RM = ^exp{-(2p 2 -l)^i}, (11) 
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Table 2: Lattice results for the slope of the Isgur-Wise function at zero recoil. 



Collab. 


Pi 


n 2 

Pu,d 


UKQCD 10 


i r, + 2 + 2 

l.Z _ 2 _ 1 


0.9 +1 ±* 


BSS 12 


1.24(26)(36) 




MO 13 


0.95 





Figure 2: UKQCD's estimate of |V c {,|£nr(w) compared to the CLEO data. 




where p 2 is the slope of at zero recoil. Lattice results for p 2 for light 
quark masses either at the strange quark mass or in the chiral limit from 
various collaborations are shown in Tab. 2. 

In order to extract |V c f,| we can now use the parametrised form of to 
extrapolate the measured decay rate for B — > D*lvt decays to u> = 1. Fitting 
the CLEO II data 14 to our estimate of |V c (,|^nr (neglecting radiative and power 
corrections), we obtain 

|V c b|= 0.037 IJ(exp) t^(stat) l^(syst). (12) 
A comparison of the experimental data and our fit is shown in Fig. 2. 



4 Heavy Baryon Spectroscopy 

Lattice studies can help making predictions for the masses of A, S, S and fl 
states in both the charm and beauty sector. Baryon masses are extracted from 
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Table 3: Heavy baryon states and their masses measured experimentally and on the lattice 16 . 



Baryon 
h = c,b 




(I) (S) si" 


Mass 
[GeV] 


Quark 
Content 


Lattice 
[GeV] 


A h 


i + 

2 


(0)(0) 0+ 


2.285(1) 
5.64(5) 


(ud)c 
(ud)b 


2.28 1 \ 
5.59 


Z h 


1 + 

2 


(1)(0) 1+ 


2.453(1) 


(uu)c 
(uu)b 


2.45 ± \ 
5.69 ±J° 




3 + 

2 


(1)(0) 1+ 


2.53(5)(?) 


(uu)c 
(uu)b 


2.43 1 4 
5.68 l x g 




1 + 

2 


(i)(-i) o + 


2.470(4) 


(us)c 
(us)b 


2.41 1 ^ 
5.69 ± X S 




1 + 

2 


(0)(-2) 1+ 


2.74(2) 


(ss)c 
(ss)b 


2.68 1 * 
5.91 1 J 



the exponential fall-off of correlation functions of suitably chosen interpolating 
operators of baryonic states. 

The results by UKQCD 15 are shown in Tab. 3 together with the experimen- 
tal numbers as far as these are known. There is good agreement between lattice 
and experiment for charmed baryons, and also the lattice value of the mass 
of the A;, is consistent. Taking the uncertainty in the lattice scale a -1 into 
account, UKQCD's result is M Ab = 5.59 t^stat) 1 gg(syst) GeV, whereas 
the PSI-CERN-Wuppertal Collaboration 16 obtain M Ab = 5.728 ± 0.144 ± 
0.018 GeV. Both these numbers are consistent with the experimental value of 
M Ab = 5.64 ± 0.05 GeV 17 . 

Mass splittings between the A, S and/or pseudoscalar mesons (P) are also 
obtained. Tab. 4 shows UKQCD's data for the dimensionless splittings Ra c b 
and i?s c6 , which show excellent agreement with the experimental numbers. 
The mass of the S* seems to be slightly below that of the S, as can be seen 
from Tab. 4. In fact, our lattice estimates for the S and S* masses are statisti- 
cally compatible. Recent measurements 18 , however, indicate that the E* — S 
splitting may be so small that it cannot be resolved with our present statistics. 

In future we will study semi-leptonic A;, — > K c lvi decays. Preliminary 
results are presented elsewhere 19 . 



6 



Table 4: UKQCD's results for the mass splittings R Ah = (M Aft - M Ph )/(M Ah + Mp h ) and 
Rs h = (-^Sh — M\ h )/(Mj; h + AfA fc ), h = c,b compared to experiment. 





Rk c R\ b 


-Rx; c -Rs c 


Lattice 
Exp. 


0.099 1? 0.033 ±1 
0.100(3) 0.033(5) 


0.038 tjj 0.017 1? 
0.035(1) 
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